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Abstract 

The predictions of event anisotropy parameters from transport model RQMD 
are compared with the recent experimental measurements for 158^4 GeV 
Pb+Pb collisions. Using the same model, we study the time evolution of 
event anisotropy at 2A GeV and 158A GeV for several colliding systems. For 
the first time, both momentum and configuration space information are stud- 
ied using the Fourier analysis of the azimuthal angular distribution. We find 
that, in the model, the initial geometry of the collision plays a dominant role 
in determining the anisotropy parameters. 



1 



Typeset using REVTgX 



Event anisotropy, often called flow, has been observed in heavy ion collisions at every 
laboratory energy [|T]-0. It is believed that the information about the equation of state 
(EOS) can be obtained from the study of flow |]IU|| . Recently, several authors argue that the 
event shape with respect to the reaction plane may carry the information about pressure 
created at the early stage of the collision [|IT|-|H)|]. 

Theoretically, in terms of hydrodynamics, the connection between the EOS and flow 
is well defined and has been studied extensively for many years. This is especially the 
case for collisions at beam energy below 15A GeV. (For the latest review of the subject, 
see Ref. |P,|16|.) However, it is still not clear whether the experimentally observed event 



anisotropy is of a dynamic origin (pressure-driven hydroexpansion) or due to shadowing 
of cold nuclear matter, passing time etc. If the initial geometry of the system plays an 
important role, the information about the collision dynamics may be obscured. In any case, 
it seems to be necessary to understand the interplay between these two competing effects. 
For this purpose, we use the transport model RQMD(v2.3 cascade mode) |DJ to study the 
event anisotropy as a function of the colliding system and beam energy. We also investigate 
the collision time dependence of the event shapes. 

In RQMD, the reaction plane is defined by the impact parameter (x direction) and 
the projectile momentum (z direction). Particle azimuthal distribution with respect to the 
reaction plane at a given rapidity window can be deconvoluted by the Fourier expansion 
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where the first and second Fourier coefficients, v\ and t>2, are connected to the directed flow 
and elliptic flow, respectively ||24|| . The coefficient vq is a normalization constant and <fi is 



defined as the azimuthal angle difference between the particle and the reaction plane. In our 
case, the reaction plane angle is explicitly set to zero. Hence angle in momentum space 
is = taxi -1 (p y /p x ), and in configuration space = taiT 1 {y / x) . At a given rapidity, the 
coefficients are: 



«i = (cos(0)>, w 2 = (cos{2(f>)). (2) 

In the following, we will first compare the recent NA49 measurements of Fourier moments 
with RQMD calculations. In section II and III, the colliding system dependence and time 
evolution of anisotropy parameters V\ and t>2 will be discussed. Finally, a brief summary 
will be presented in section IV. 

I. COMPARISON OF DATA WITH MODEL PREDICTIONS FOR 1584 GEV 

PB+PB COLLISIONS 

The experimentally measured (fi,^) parameters for proton and pion as a function of 
rapidity are shown as filled circles in Figure |l|. These data are extracted from the 1584 GeV 
Pb+Pb collisions by the NA49 collaboration |||7]]. More details about methods involved in 



the experimental determination of these parameters can be found in the recent paper |19 
The model calculations (open circles) are done within similar impact parameter range and 
acceptance (pt,y) cuts. As one can see, the model calculations are in reasonable agreement 
with the data for both protons and pions. The sign of the first Fourier moment v\ for pions 
is opposite to that of protons for both data and model calculations indicating the nuclear 
shadowing effects in the collision. An interesting feature in Figure [l] is that there seems to be 
two slopes for the function Vi(y): at \y\ < 1.5, the distribution is relatively flatter compared 
to the region \y\ > 1.5. This implies that the physics around the mid-rapidity region is 
different from that near the projectile (target) rapidity region. One possible interpretation 
is that near the beam rapidity region, the number of rescatterings is small, therefore the 
initial geometry leads to a large value of V\. On the contrary, a large number of rescatterings 
washes out the initial memory and leaves small anisotropy around the mid-rapidity. 



As predicted in Ref. ||20|| , the second Fourier moment t>2 is found to be positive for both 
nucleons and pions at this beam energy. This implies an enhancement of in-plane emission 
from the particle source. More discussions on this point will be presented later in this paper. 

Since the model predictions are in reasonable agreement with the experimental data, it 



may be instructive to study system size dependence and the space-time evolution of the 
anisotropy parameters within the model. 



II. Vi FOR DIFFERENT COLLISION SYSTEMS 

In Figure 0, the model predicted v\ for nucleons and pions as a function of rapidity are 
shown for 2 A GeV S+S (3 < b < 5 fm), 2 A GeV Ru+Ru Q (4 < b < 6 fm), 2 A GeV Pb+Pb 
(5 < b < 8 fm), 158 A GeV S+S (3 < b < 5 fm) and 158A GeV Pb+Pb (5 < b < 8 fm) 
collisions. The impact parameter range is chosen such that the maximum value of V\ is 
achieved for any given colliding system. 

It is known from the hydrodynamic calculations that the pressure created at the early 
stage will lead to sizable values of V\ and t>2, and the pressure is different for different bom- 



barding energies |[21|| . However, from Figure |2| and Table |, one notices that the maximum 
values of v i for nucleons are independent of the beam energy for a given colliding system. 
For a fixed beam energy, the maximum value of V\ depends on the size of the colliding nuclei. 
The larger the size, the higher the maximum value of V\. This may suggest that the initial 
geometry, i.e. the size of the colliding nuclei and impact parameter, determines the strength 
of V\ near the beam (target) rapidity region, providing the elastic scattering amplitude does 
not vary much within energy range discussed here. Similar trends are also seen for pions 
even though the V\ for pions has the opposite sign to that of nucleons due to the shadowing 
effect. 

The other interesting feature is that the slope of the v i curve around the mid-rapidity 
region depends on beam energy for both nucleons and pions. The higher the beam energy, 
the smaller the slopeQ. The value of v 2, on the other hand, is maximized around mid-rapidity. 
It was suggested that the second Fourier moment may carry the information about the early 



x The exact beam energy for the Ru+Ru collisions is 1.69^4 GeV. 
2 This is also true for v\ as a function of normalized rapidity. 
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stage of the collision [11], 15]. The sensitivity of the V2 signal to the collision dynamics will 
be discussed in the next section. 



III. TIME EVOLUTION OF (Vi,V 2 ) 

Now let us explore (vi, v 2) of nucleons and pions as a function of freeze-out time for two 
colliding systems: 2A GeV Au+Au and 158A GeV Pb+Pb collisions. In a discussion of the 
time evolution, a brief remark concerning the nuclear passing time is necessary. Nuclear 
passing time is the time needed for two incoming nuclei to pass each other and defined as 
t p = 2r/(7/3) where r is the nuclear radius. Passing time sets the scale on which the effects of 
possible nuclear shadowing, e.g. the interaction of the participating nucleons and produced 
particles with the spectator matter, are significant. Obviously, the passing time depends 
on the beam energy and on the size of colliding nuclei. The nuclear passing time t p , in the 
nucleon-nucleon center-of-mass frame, for 2A GeV Au+Au and 158A GeV Pb+Pb collisions 
are listed in Table 0. Also listed are the nucleon mean freeze-out time t m from the RQMD 
calculations. The ratio of t m /t p for the low energy collisions is much smaller than that at 
the high energy, indicating that the nuclear shadowing is relatively more important for lower 
energy collisions. 

A. vi,V2 (momentum space) 

The calculated (t+t^) values, determined in the momentum space, as a function of 
rapidity in different time windows are shown in Figures [3] and f| for nucleons and pions, 
respectively. In these figures, the open circles are for the lower energy (2 A GeV) and the 
filled circles for the higher energy (158^4 GeV) collisions. Similar distributions calculated at 
the corresponding passing time (0 < t < 12 fm/c for 2A GeV and < t < 2 fm/c for 158A 
GeV collisions) are shown in Figure [5] 

For nucleons, as one can see in Figure |^ (left column lower 4 plots), the strength of v± 
is gradually built up as a function of collision time. At the earlier time, t < 10 fm/c, a 
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well defined negative slope (around mid-rapidity) is seen for the low energy collisions while 
at the high energy the distribution is relatively flat. This flatness is due to the averaging 
over the 10 fm/c time period which is much longer than the passing time for the 158A GeV 
collisions (see Table [TT]). In fact, the negative slope is also present at that energy for t < t p 
as shown in Figure This again suggests that the shadowing of the cold nuclear matter 
is responsible for the observed negative slope at the early stage of the collision. At t > 10 
fm/c, the maximum value of v\ has been well established for the lower energy collisions. For 
the higher energy collisions, the fully established v\ only appears after t > 30 fm/c. 

The time dependence of v\ for pions shows a different trend (Figure 4 left column). The 
maximum values of v i are achieved at the early stage of the collision. As time goes on, the 
distributions of v % become flatter. The fact that the v\ distributions for pions always have 
a negative slope through the whole collision history suggests the importance of shadowing 
effect. However, the shadowing matter is different for pions at different collision stages. At 
early time, t < t p , pions (as well as the participating nucleons) are shadowed by the cold 
spectators. Later, after the spectator matter leaves the collision zone, pions are shadowed by 
the participant nucleons. This may be the underlying mechanism that leads to the different 
behavior of v\ for nucleons and pions. 

The second Fourier coefficient describes an elliptical shape of the azimuthal distribu- 
tion. A positive sign of v-i means that the longer axis of the ellipse is oriented in the reaction 
plane while a negative sign of v% indicates the longer axis is in the direction perpendicular 
to the reaction plane. This latter effect was called "squeeze-out" at Bevalac energies p2| . 



In Figure ^ and Figure ^ (right column), the t>2 distributions are shown as a function of 
rapidity in different time intervals for nucleons and pions, respectively. It can be seen that 
the time averaged values of V2 for nucleons and pions at both beam energies are positive, 
implying that particles preferentially freeze-out in the reaction plane. 

However, at the time scale comparable with nuclear passing time, nucleons and pions 
demonstrate an out-of-plane emission pattern ("squeeze-out") for 2 A GeV collisions (Fig- 
ure fright columns) which is consistent with the cold nuclear matter shadowing scenario. For 
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high energy collisions, the "squeeze-out" effect is absent for both pions and nucleons. This 
may be explained by a much shorter nuclear passing time and higher Lorentz contraction 
of the colliding nuclei. According to the RQMD calculations, at about the corresponding 
passing time, the average number of collisions suffered by the mid-rapidity nucleons is 3 
and 5 for 158 A GeV and 2 A GeV collisions, respectively. In the high energy collisions the 
rapidity span is about 6 units. Hence, three collisions are just about enough to move a 
nucleon from projectile rapidity to mid-rapidity. At this energy, neither time nor number of 
rescatterings is sufficient to develop an out-of-plane emission. 



B. t>i,t>2 (configuration space) 

Up to now, we have discussed the anisotropy parameters (v i, t> 2 ) in the momentum space. 
It is interesting to look at these parameters in coordinate space in order to understand the 
space-time evolution of the collisions within the model. To the best of our knowledge such 
studies have not been performed before. The information about source anisotropy may be 



obtained experimentally by studying two-particle correlation functions |23|-p6 . 

In general, the v\ parameters in the configuration space for both nucleons and pions (see 
Figure |6] and |7| left columns) behave differently from those in the momentum space: (1) 
The shapes of the v i distributions are different and the maximum values of V\ (around the 
beam and target rapidities) in the configuration space are much larger than those in the 
momentum space; (2) The v\ of pions has the same sign to that of nucleons; (3) The v\ 
becomes stronger as a function of time, especially for pions. Furthermore, similar to what 
have been observed in the momentum space, the maximum values of V\ are the same for 
the two collisions at different beam energies. This implies that the initial collision geometry 
plays an important role here. 

For the non-zero impact parameter heavy-ion collisions, the transverse shape of the 
initial collision zone is an ellipse with its longer axis oriented perpendicular to the reaction 
plane (t> 2 < 0). Such an event shape is indeed seen for both collision energies in Figure 
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(right column) for the time interval t < 10 fm/c. As time goes on, rescattering occurs among 
particles which leads to a more spherical shape of the reaction zone, v 2 — > 0. Note one would 
expect a large positive v 2 at later times in the case of pressure driven hydroexpansion. For 
the time integrated distributions, some remnants of the initial shape still can be seen for the 
2A GeV collisions while for the higher energy collision the source appears to be isotropic. 

IV. SUMMARY 

Based on the RQMD(v2.3) model calculations, we studied the anisotropy parameters 
(vi,v 2 ) as a function of time for different energies and different colliding systems. From a 
comparison between the model calculations and the recent experiment measurement at the 
beam energy of 158A GeV, we find that the model provides a reasonable description of the 
observed event anisotropy. From the study of energy dependence, we find that the combi- 
nation of collision geometry and nuclear passing time plays a dominant role in determining 
the behavior of the v± and v 2 parameters. Our studies show that the anisotropy parameters 
extracted in the momentum space are different from the ones extracted in the configuration 
space. As a function of collision time, some rather dramatic changes are found for the event 
anisotropy parameters in both momentum and configuration spaces. At t < t p , the pattern 
of the distributions is strongly influenced by the spectators. At later times, a prominent 
in-plane emission is evident. 
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TABLES 



Colliding System 


A 


r (fm) 


y{ V ™ ax ) 


„,max 
v l 


2A GeV S+S 


32 


3.8 


0.91 


0.171+0.007 


2A GeV Ru+Ru 


96 


5.5 


0.91 


0.186+0.003 


2A GeV Pb+Pb 


208 


7.1 


0.91 


0.225+0.009 


158A GeV S+S 


32 


3.8 


2.92 


0.166+0.008 


158A GeV Pb+Pb 


208 


7.1 


2.92 


0.225+0.009 



TABLE I. Nucleon v\ parameters for symmetric A — A collisions, the atomic number A, the 
radius (r = 1.2 • A 1 / 3 ) of the beam(target) nucleus and the rapidity corresponding to the maximum 
value of v\ are shown. 



Energy 


/3(c) 


7 


y 


t p (fm/c) 


t m (fm/c) 


2 GeV 


0.718 


1.438 


0.904 


13.52 


20.63 


158 GeV 


0.994 


9.231 


2.913 


1.55 


33.89 



TABLE II. Velocity (3 (in the unit of speed of light), 7, rapidity y and passing time t p are 
shown for 2A GeV Au+Au and 158A GeV Pb+Pb collisions. The RQMD mean freeze- out time t m 
for nucleons are also included. All quantities are evaluated in the nucleon-nucleon center- of-mass 
frame. 
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FIG. 1. v\,V2 for proton and pion as a function of rapidity (filled circles - experimental results, 
open circles - model calculations) from 158^4 GeV Pb+Pb collisions. 
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FIG. 2. v\ as a function of rapidity y, in center-of-mass frame, for nucleons and pions for 
different colliding systems. The arrows indicate the beam and target rapidities. 
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FIG. 3. The time evolution of v\ and V2 for nucleons in the momentum space as a function of 
rapidity in the center-of-mass frame for 2A GeV Au+Au collisions (open circles) and 158A GeV 
Pb+Pb collisions (filled circles). 
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FIG. 4. The time evolution of v\ and V2 for pions in the momentum space as a function of 
rapidity in the center-of-mass frame for 2A GeV Au+Au collisions (open circles) and 158A GeV 
Pb+Pb collisions (filled circles). 
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FIG. 5. Nucleon and pion anisotropic parameters V\,V2 as a function of rapidity, around the 
nuclear passing time, for 2A GeV Au+Au (open circles) and 158.A GeV Pb+Pb (filled circles) 
collisions. 
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FIG. 6. The time evolution of v\ and V2 for nucleons in the configuration space as a function 
of rapidity in the center-of-mass frame for 2A GeV Au+Au collisions (open circles) and 158A GeV 
Pb+Pb collisions (filled circles). 
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FIG. 7. The time evolution of v\ and V2 for pions in the configuration space as a function of 
rapidity in the center-of-mass frame for 2A GeV Au+Au collisions (open circles) and 158A GeV 
Pb+Pb collisions (filled circles). 
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